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ABSTRACT 

Aims. This is the first paper of a series aimed to understand the formation and evolution of bars in early-type spirals and their 
influence in the evolution of the galaxy. 

Methods. Optical long-slit spectra along the major axis of the bar of a sample of 20 galaxies are analyzed, velocity and velocity 
dispersion profiles along the bar are presented. Line-strength indices in the bar region are also measured to derive stellar mean- 
age and metallicity distributions along the bars using stellar population models. 

Results. We obtain mean ages, metallicities and chemical abundances along the bar of 20 galaxies with morphological types from 
SB0 to SBbc. The main result is that we find large variation in age and metallicity along the bar in 45% of our sample. We find 
three different types of bars according to their metallicity and age distribution along the radius: 1) Bars with negative metallicity 
gradients. They show mean young/intermediate population (< 2 Gyr), and have amongst the lowest stellar maximum central 
velocity dispersion of the sample. 2) Bars with null metallicity gradients. These galaxies that do not show any gradient in their 
metallicity distribution along the bar and have negative age gradients (i.e younger populations at the bar end). 3) Bars with 
positive metallicity gradients, i.e. more metal rich at the bar ends. These galaxies are predominantly those with higher velocity 
dispersion and older mean population. We found no significant correlation between the age and metallicity distribution, and 
bar/galaxy parameters such as the AGN presence, size or the bar strength. From the kinematics, we find that all the galaxies 
show a disk-like central component. 

Conclusions. The results from the metallicity and age gradients indicate that most galaxies with high central stellar velocity 
dispersion host bars that could have been formed more than 3 Gyrs ago, while galaxies with lower central velocity dispersions 
show a wider distribution in their population and age gradients. A few bars show characteristics compatible with having been 
formed less than <2 Gy ago. However, we do not have a definite answer to explain the observed gradients and these results 
place strong constrains to models of bar formation and evolution. The distribution of mean stellar population parameters in 
the bar with respect to a is similar to that found in bulges, indicating a close link in the evolution of both components. The 
disk-like central components also show the important role played by bars in the secular evolution of the central structure. 

Key words. Galaxies: spirals- Galaxies: abundances - Galaxies: evolution - Galaxies: formation - Galaxies: structure - Galaxies: 
kinematics 
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1. Introduction 

Bars are known to be efficient mechanisms to redistribute 
angular momentum and matter in the galaxies up to large 
distances and, therefore, they are expected to play an im- 
portant role in the secular evolution of disk galaxies. The 
bar gravity torques make the gas lose angular momen- 
tum which, in turn, provokes an inflow of gas towards 
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the central parts. This inflow of gas accumulates in cen- 
tral mass concentrations (CMC), fueling the central black 
hole (e.g Sholsman et al. 1989) and, possibly, forming a 
stellar bulge. Bars are found in ~60% of spiral galaxies 
in the local Universe (e.g. Eskridge et al. 2000; Whyte 
et al. 2002; Menendez-Delmestre et al. 2007; Marinova & 
Jogee 2007), or 30% if only strong bars are considered. 
Despite their obvious importance in disk galaxy evolution, 
their own fate is still a matter of debate. Several mecha- 
nisms are able to destroy the bar: for example; 1) the bar 
may dissolve in the presence of a sufficiently massive cen- 
tral component (e.g. a black hole) (Friedli & Pfenniger 
1991; Friedli & Benz 1993, Norman, Sellwood & Hasan 
1996; Berentzen et al. 1998; Fukuda et al. 2000; Bornaud 
& Combes 2002; Shen & Sellwood 2004; Athanassoula, 
Lambert & Dehnen 2005; Bornaud, Combes & Semelin 
2005). However, the masses of the CMC needed to de- 
stroy the bar are much larger than those of the most su- 
per massive black holes found so far in disk galaxies (Shen 
& Selwood 2004; Athanassoula, Lambert & Dehnen 2005; 
Ferrarese & Ford 2005). 2) Transfer of angular momen- 
tum from infalling gas to the bar has also been proposed 
as a mechanism that could strongly weaken the bar. The 
combined effect of this angular momentum transfer and 
the CMC can destroy the bar in Sb-Sc galaxies in 1-2 Gyr 
(Bournaud & Combes 2002). Then, if gas is present, a new 
bar could form with characteristics different from those of 
the parent bar. However, simulations of the effect of gas 
on the stellar evolution of a bar embedded in a live halo 
(Berentzen et al. 2007) have shown that, althogh there are 
some structural differences in bars evolving with and with- 
out gas, in both cases the stellar bar can survives, at least, 
5 Gyr (the whole computation time of these simulations 

Do all bars go trough these processes? Do some bars die 
while others are robust over many Hubble times? Previous 
studies have tried to put some limits to the age of the 
bars by analyzing the evolution with redshift of the bar's 
fraction in galaxies. These studies have also obtained dis- 
crepant results: Some authors have found a fraction of 
barred galaxies at z>0.5 lower than the local fraction 
(Abraham et al. 1999), although some authors claim that 
this may be the consequence of selection effects, due to the 
high angular resolution needed to find bars (Elmegreen, 
Elmegreen & Hirst. 2004; although see van den Bergh et al. 
2002). In order to avoid the problem of the angular reso- 
lution, several studies have carried out this analysis using 
the ACS camera on the HST: Elmegreen et al. (2004) and 
Jogee et al. (2004) found the same bar fraction (~0.4) at 
redshift z=l.l than in the local Universe, suggesting that 
the bar dissolution cannot be common during a Hubble 
time unless of course, the bar formation rate is compa- 
rable to the bar destruction rate. On the contrary, Sheth 
et al. (2008), in a recent study, using images from the 
COSMOS survey (Scoville et al. 2007) and using a larger 
sample than previous studies, found that the bar fraction 
at z=0.84 is one-third of the present-day value. They also 
found a much stronger evolution for low mass galaxies and 
late-type spectral types. Part of the discrepant results may 



be due to the selection effects and other systematic effects 
that need to be still further investigated. 

Detailed analysis of the stellar populations in the bar 
region of local galaxies can shed some light to the for- 
mation and evolution of bars. Gadotti & de Souza (2006) 
obtained the color and color gradients in the bar region of 
a sample of 18 barred galaxies. They interpreted the color 
differences as differences in stellar ages. They concluded 
that younger bars were hosted by galaxies of later types. 
However, the conclusions are hampered by the assumption 
that the color traces the age, and that the metallicity ef- 
fect or dust extinction are not important. Clearly, a study 
attempting to break the age-metallicity degeneracy and 
taking into account the effect of dust is needed. In ad- 
dition, the influence of dust in the line-strenght indices 
have been shown to be minimum (MacArthur 2005) while 
ages and metallicities from broad-band colors are heavily 
affected by dust as well as by the age-metallicity degen- 
eracy. 

To date, it has been difficult to obtain stellar popu- 
lation parameters along the bar because a high signal-to- 
noise is required for this analysis. Although bars are high 
surface brightness structures, this still implies long inte- 
gration times on medium-size telescopes. Recent work has 
presented the radial distribution of line-strength indices 
along bars (Perez, Sanchez-Blazquez & Zurita 2007) for 6 
galaxies. 

The theoretical model of abundance gradients in bars 
developed by Friedli (1994) predicted, using N-body simu- 
lations of bars with pre-existing exponential abundances, 
a null evolution of the stellar abundance profile along the 
bar (although with a decrease in the mean metallicity) 
while the gas abundance profile flattened rapidly. Studies 
of the radial gas-phase abundance distribution have shown 
that there is, indeed, a flattening in the gas abundance 
along the bar (Martin & Friedli 1997, 1999). However, re- 
cent work (Perez, Sanchez-Blazquez & Zurita 2007) has 
shown, from stellar absorption line-strengths, that some 
bars present stellar metallicity gradients opposite to those 
found in the disks (i.e. more metal rich population at the 
end of the bar than in the regions of the bar closer to the 
center) , contradicting the numerical results. This surpris- 
ing result needed a more systematic study of the stellar 
population along the bars, properly comparing the line- 
strength indices of a larger sample of galaxies with the 
state-of-the art stellar populaton models to derive and 
quantify the variations of age and metallicity along the 
bar. 

In this paper, we present a detailed stellar population 
analysis of the bar region of a sample of 20 early-type spi- 
rals (from SO to Sbc) with and without nuclear activity 
in order to derive ages, metallicities and relative chemical 
abundances ratios. The sample selection is described in 
Sect. [5] The observations and data reduction are presented 
in Sect. [3l The velocity dispersion measurements and the 
emission correction applied to the data are explained in 
Section 14.11 The line-strength derivation is described in 
Sect. 14.21 The morphological characterization of the sam- 
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pie is described in Sect. 14.31 Section |4~41 gives a description 
of the methodology followed to derived the stellar popu- 
lation parameters. The results of the line strength-indices 
and the kinematics are presented in Sect. In this work, 
we concentrate only in the stellar populations of the bar 
region. The results of the central indices, bulge and bulge 
gradients will be presented in the second series of this pa- 
per. 

2. Sample characterization 

We have selected barred galaxies from the Third Reference 
Catalogue of bright galaxies (RC3), de Vaucouleurs 1948, 
with the following criteria; to be classified as barred 
and with inclinations between 10° and 70°, and nearby 
(cz < 4000 kms -1 ) to be able to properly resolve the bar. 
The sample is biased towards early-types barred galax- 
ies, which have higher surface brightness, to ensure a high 
signal-to-noise ratio (S/N), crucial for a reliable determi- 
nation of line-strength indices. Since we also aim to relate 
the mean stellar ages and metallicities of the bar to the 
nuclear activity in the center of the galaxy, half of the 
galaxies were selected to have nuclear activity. Eight of 
the galaxies present nuclear bars. Our final sample com- 
prises 20 galaxies. It is, by no means, a statistically com- 
plete sample, but it is a representative sample of early-type 
galaxies bar population. Table[T]shows the main properties 
of the sample as taken from the Hyperleda galaxy cata- 
log (Paturel et al. 2003) Q. The bar strengths (shown in 
Tab.Q]as bar-class; for assigment of bar-class to a certain 
bar strength, see Buta & Block 2001, is a scale from to 6, 
being 6 the strongest bar) have been taken from the liter- 
ature, where strength is defined as the torque of a bar em- 
bedded in its disk (Combes & Sanders 1981), see TableQ] 
for the references for the individ ual galaxies. The nu- 
clear t ypes have been obtained from lVeron-Cettv fc Veronl 
( 2006f ). The sample shows a wide distribution of maximum 
rotational velocities (80 - 260 kms -1 ). 

3. Observations and data reduction 

We obtained long-slit spectra along the bar in our sample 
of 20 barred galaxies. The bar position angles were derived 
using the Digital Sky Survey (DSS) images. 

The observations were performed in two different runs. 
In the first run (runl hereafter), spectroscopy of 6 galax- 
ies was obtained at the 2.3m telescope at Siding Spring 
Observatory (RSAA, ANU, Australia) during February 
2006 with the Double Beam Spectrograph (DBS, Rodgers, 
Conroy & Bloxham 1998). The gratings 600B and the 
600R were used for the blue and red arms, respectively, 
with a slit-width of 2 arcsec. This set-up gives a disper- 
sion of 1.1 A pixel -1 for the blue-arm and 1.09 A pixel -1 
for the red arm in the wavelength intervals 3892-5814 A 
and 5390-7314 A respectively, and a final spectral resolu- 
tion of FWHM- 2.2 A. One of the main problems in the 



http:/ /leda.univ-lyonl.fr 



derivation of accurate line-strength indices in regions of 
low surface brightness comes from the systematic effects 
due to sky subtraction. One of the advantages of the DBS 
is the length of the slit (6.7 arcmin) which allows the se- 
lection of sky regions at large radii, not contaminated by 
the light of the galaxy. 

The second set of observations (run2 hereafter) were 
carried out in February 2007 at El Roque de Los 
Muchachos Observatory (La Palma, Spain). We used the 
IDS mounted on the Isaac Newton telescope (2.5m). The 
R632V grating combined with a 1.5 arcsec slit-width gives 
a resolution of 3 A (FWHM). The full unvignetted slit 
length in IDS is 3.3 arcmin. In order to obtain enough 
sky coverage in one of the galaxy sides, we displaced the 
center of the slit from the galaxy center during the obser- 
vations. 

Comparison arc lamp exposures were obtained 
in both runs to wavelength calibrate the frames. 
Spectrophotometric standards were observed with a slit- 
width of 6 arcsec in order to avoid differential flux-loses 
due to atmospheric diffraction. Additionally, we observed 
11 and 35 G-K stars in the first and second run respec- 
tively from the Lick/IDS (Gorgas et al. 1993; Worthey 
et al. 1994) and MILES libraries (Sanchez-Blazquez et al. 
2006) in order to calibrate our data into the Lick/IDS 
spectrophotometric system. The total integration time for 
each galaxy is listed in Table [51 Only the blue part of the 
spectra in the case of the DBS observations is analyzed in 
this paper. The red part of the spectra, i.e. the nebular 
emission properties, will be presented elsewhere. 

The reduction of the two runs has been carried out 
with the package REDUCEME (Cardiel 1999). Standard 
data reduction procedures (flat-fielding, cosmic ray re- 
moval, wavelength calibration, sky subtraction and flux- 
ing) were performed. Error images were created at the 
beginning of the reduction and were processed in paral- 
lel with the science images. Initial reduction of the CCD 
frames involved bias subtraction and removal of pixel- 
to-pixel sensitivity variations (using flat-field exposures 
of a tungsten calibration lamp). Correction of a two- 
dimensional low-frequency scale sensitivity variations was 
done using twilight sky exposures. Arc frames were used 
to convert the spectra into a linear wavelength scale. We 
use typically ~ 60 lines. The lines were fitted with a 1st 
order polynomia, with r.m.s. errors lower than 0.1 A. 
Atmospheric extinction was calculated using the extinc- 
tion curve of King (1985). To correct the effect of inter- 
stellar extinction, we used the curve of Cardelli, Clayton 
& Mathis (1989). The reddenings applie d were extracted 
from the RC3 catalogue of galaxies ( de Vaucouleursl 
Il948 ). Relative flux calibration was achieved using the 
spectra of the observed standard stars. All the flux- 
calibration curves of each night were averaged and the 
flux calibration errors were estimated by the differences 
between the indices measured with different curves. After 
correction for geometrical distortions (both spatial and 
spectral) , a sky image was generated for each observation 
by fitting, for each channel, a low order polynomial using 
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Table 1. General properties of the sample 



Object 


v(km s i ) 


Type 


Bar-class 1 


Nuclear type 


Inner morph.' ! 


B 


^max,gas (km S ) 


i(deg) 


NGC 1169 a 


2387 


SABb 


3 






12.35 


259.1 ± 7.3 


57.1 


NGC 1358 


4028 


SAB(R)0 


— 


Sy2 


— 


13.19 


136.1 ± 10.6 


62.8 


NGC 1433 6 


1075 


(R)SB(rs)ab 


4 


Sy2 


Double-bar" 


10.81 


85.1 ± 2.4 


68.1 


NGC 1530 6 


2461 


SBb 


6 


_ 


_ 


12.50 


169.1±3.5 


58.3 


NGC 1832 d 


1939 


SB(r)bc 


2 


_ 


_ 


12.50 


129.9±2.0 


71.8 


NGC 2217 


1619 


(R)SB(rs)0/a 


_ 


LINER? 


Double-bar 6 


11.36 


183.4± 9.2 


30.7 


NGC 2273 c 


1840 


SB(r)a 


2 


Sy2 


_ 


12.62 


192.2±5.5 


57.3 


NGC 2523 


3471 


SBbc 






_ 


12.64 


211.4 ±10.9 


61.3 


NGC 2665 


1734 


(R)SB(r)a 


_ 


_ 


_ 


12.47 


130.9±7.1 


32.8 


NGC 2681 c 


692 


(R)SAB(rs)0/a 


1 


Sy3 


Triple-bar c 


11.15 


87.5±6.7 


15.9 


NGC 2859 c 


1687 


(R)SB(r)0 " 


1 


Sy 


Double-bar d 


11.86 


238.5±13.3 


33.0 


NGC 2935 


2271 


(R)SAB(s)b 


_ 




_ 


12.26 


188.3±2.0 


42.7 


NGC 2950 


1337 


(R)SB(r)0" 


_ 





Double-bar e 


11.93 


_ 


62.0 


NGC 2962 


1966 


(R)SAB(rs)0 


- 


- 


Double-bar c 


12.91 


202.9±9.9 


72.7 


NGC 3081 c 


2391 


(R)SAB(r)0/a 


3 


Sy2 


Double-bar d 


12.89 


99.9±4.0 


60.1 


NGC 4245 c 


815 


SB(r)0/a 


2 






12.33 


113.5±5.4 


56.1 


NGC 4314 a 


963 


SB(rs)a 


3 


LINER 


Double-bar c 


11.42 


253.3±24.6 


16.2 


NGC 4394 d 


922 


(R)SB(r)b 


3 


LINER 




11.59 


212.5±16.0 


20.0 


NGC 4643 c 


1335 


SB(rs)0/a 


3 


LINER 




11.68 


171.4±7.2 


42.9 


NGC 5101 d 


1868 


(R)SB(r)0/a 


2 






11.59 


195.7±9.0 


23.2 



(1) : 

(a) Bar class derived from the K — band light distribution. 
(6) Bar class derived from the K — band light distribution, 

(c) Bar class derived from the K — band light distribution, 

(d) Bar class derived fro m the H — band l ight d istri bution , 

(2) : (a) iButal (|l986T ) (6) lJungwiert et al.l (1997T ) ( c) lErwinl 

(3) Rotational velocity corrected for inclination 



iBlock et all d200ll ) 
Block et al.r(l2004' 



Buta et al 



Laurikaincn et al. (2004)) 

(|2004l ) (d) IWozniak et al l (|l995l ) 



regions selected at both sides of the galaxy frame in the 
first run, and only in one of the sides of the frame for the 
galaxies in the second run. 

In order to derive the index spatial distribution, for 
the fully reduced galaxy frames, a final frame was created 
by extracting spectra along the slit, binning in the spatial 
direction to guarantee a minimum signal-to-noise ratio of 
20 per A in the spectral region of Mgb. This minimum 
S/N ensures errors lower than 15% in most of the Lick/IDS 
indices (Cardiel et al. 1998). 

4. Analysis 

4.1. Velocity dispersion measurements and emission 
correction 

Whenever ionized gas is present in the galaxy, emission 
lines can affect the measurement of the absorption in- 
dices. In particular, all the Balmer lines (H.8, H7 and H/3), 
Fe5015, and Mgb indices are affected by the presence of 
H5A4101 , H 7 A4340, H/3A4861, [OIII]AA 4959,5007, and 
[NII]AA5198,5200 lines in emission. The presence of emis- 
sion in the Balmer lines tend to fill the absorption pro- 
files making the indices weaker. If this is not corrected 
properly, older ages can be artificially measured. In order 
to correct from the emission lines contamination we have 
used the routine GANDALF (Sarzi et al. 2005, S05 here- 



after). This routine fits, simultaneously, the stellar and the 
emission line spectra, by treating the emission lines as ad- 
ditional Gaussian templates and iteratively searching for 
the best radial velocity and velocity dispersion. To model 
the stellar spectra, we use the stellar population models 
by Vazdekis et al. (2008) based on the MILES stellar li- 
brary (Sanchez-Blazquez et al. 2006), previously degraded 
to our instrumental resolution. The continuum shape of 
the stellar templates and that of the galaxies is fitted and 
removed using a multiplicative Legrende polynomial and 
subtracted from all spectra. 

As it is explained in S05, despite that the kinematic 
of the different emission lines can be fitted independently, 
in some cases, as in H/3, degeneracies between the stel- 
lar and emission templates can cause spurious detections. 
Furthermore, it is very difficult to constrain the kinematic 
of the weak emission lines. For these reasons, we follow 
the same procedure as in S05 and impose that all the 
emission lines have the same kinematic as the strong line 
[OIII]A5007, except for the cases when H/3 line is very 
strong. The steps followed are: (1) Masking the emission 
lines and calculating a first estimate of the stellar veloc- 
ity and velocity dispersion. (2) Then, convolving our tem- 
plates with the calculated velocity dispersion and shifting 
them according to the recessional velocity, and then, cal- 
culating the kinematics of the [OIII]A5007 lines. (3) We 
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force all the emission lines to have the same kinematic 
as the one calculated with [OIIIJA5007 and fit, simultane- 
ously, emission and absorption spectra. 

Figure [1] shows the central spectrum of NGC1358, a 
galaxy with particularly strong emission, before and after 
cleaning the emission using the method described above. 
In some cases (i.e in the nuclei of active galaxies) , the emis- 
sion lines cannot be fitted with a single Gaussian, as the 
emission line profiles show asymmetries. In those cases, 
double Gaussians are fitted in order to reproduce the ob- 
served profile. This process will be explained in more detail 
in a future paper, when we analyze the properties of the 
bulges and central parts in our sample (Perez, Sanchez- 
Blazquez, & Zurita 2009, in preparation). 

GANDALF makes use of PPXF (penalized pixel- 
fitting) to calculate the stellar kinematics. PPXF is a soft- 
ware described in Capellari & Emsellem (2004). It is spe- 
cially appropriate for low signal-to-noise data or when the 
data is not well sampled. The errors in the radial velocity 
and velocity dispersion were computed through numeri- 
cal simulations. In each simulation, a bootstrapped galaxy 
spectrum, obtained using the error spectrum provided by 
the reduction with REDUCEME, is fed into the algorithm. 
Errors in the parameters are then calculated as the unbi- 
ased deviation of the different solutions. The kinematic of 
the gas and the diagnostic emission lines will be subject 
of a future paper (Florido et al. 2009, in preparation). 

4.2. Line-strength indices 

Lick/IDS line-strength indices using the definition in 
Trager et al. (1998) are measured in all the -cleaned from 
emission lines- binned spectra along the radius. The er- 
rors are calculated by the uncertainties caused by photon 
noise, wavelength calibration and flux calibration. 

Line-strength indices depend on the broadening of the 
lines caused by instrumental resolution and by the inter- 
nal motions of the stars. Therefore, when comparing with 
the stellar population models, the synthetic and observed 
indices need to be measured at the same resolution. 

We are using two different sets of models to derive 
the synthetic indices. The first one, Vazdekis et al. (2008) 
is built using a library (MILES) with constant resolution 
along the wavelength direction (2.3 A) which is also flux 
calibrated (Sanchez-Blazquez et al. 2006). The main ad- 
vantage of these models is that they predict the whole 
spectrum for a population of a given age and metallic- 
ity and, therefore, it can be convolved with a broadening 
function to mimic the desired resolution. To compare with 
these models we degrade all the spectra to the resolution 
of the runl spectra (FWHM«5.2 A). 

The second set of models we are using are those of 
Thomas, Maraston & Bender (2003, TMB03 hereafter). 
These models are based on the Lick/IDS fitting functions 
(Gorgas et al. 1993; Worthey et al. 1994) and, therefore, 
their resolution cannot be modified. To compare our mea- 
surements with this model we degrade our spectra to the 



wavelength-dependent resolution of the Lick stars follow- 
ing the prescriptions given in Worthey & Ottaviani (1997). 

Absorption lines are not only broadened by the line- 
spread function of the instrument, but also by Doppler 
broadening due to the velocity dispersion of the stars. 
To compare with the models and with galaxies with dif- 
ferent velocity dispersions, this needs to be corrected. 
Because the variation of the indices with the velocity dis- 
persion depends on the strength of the index itself (e.g. 
Kuntschncr 2000), one should perform a different correc- 
tion for each individual spectrum along the radius. To 
calculate these corrections we use the template obtained 
in GANDALF when fitting the spectra in the first bin, i.e. 
the bar-bin closest to the galaxy center. We degrade then 
this spectrum at the Lick/IDS resolution and broaden 
it to different velocity dispersions in steps of 10 kms^ 1 . 
We then build polynomia (for each individual spectrum) 
of the form: C(a)=l(a=0) /1(a) for atomic indices, and 
C(o-)=I(er = O)-I(cr) for molecular indice, where I(<r=0) 
represents the index at the Lick resolution and I(cr), the 
index at Lick resolution broadened with a velocity disper- 
sion a. Polynomia arc built as well with the initial broad- 
ening equal to the instrumental resolution of runl, in order 
to compare with the V08 models at the resolution of the 
data. We finally apply a single correction to each galaxy 
using the polynomia derived for the spectra at the begin- 
ning of the bar, and use the rest of the templates along 
the bar to calculate the error for this approximation. This 
error is added in quadrature to the final index error. 

Once the indices in the galaxy spectra and models 
are measured at the same resolution and corrected by 
the velocity dispersion broadening, a further correction is 
needed in order to fully transform the measurements into 
the right spectrophotometric system. These small correc- 
tions try to compensate for the differences between the 
models and the data that appear due to a deficient cal- 
ibration of the Lick/IDS stellar spectra. To obtain these 
offsets, we compare the index measurements in our sample 
of stars in common with the Lick/IDS library observed for 
this purpose. The final offsets are very small and listed in 
Table IB. II Figures with the comparison of individual in- 
dices can be found in appendix|Bj In principle, because our 
data are relatively flux calibrated, a further offset is not 
necessary to transform our measurements into the MILES 
spectrophotometric system (that is also flux calibrated). 
However, we also compare the indices measured in stars 
in common with this library. As expected, for most of the 
indices the offsets are null (see Table [B~Tj) . although small 
offsets were found for some of them. 

4.3. Morphological characterization 

In order to relate the stellar population properties with 
the morphological features, we perform ellipse fitting anal- 
ysis on available archived R — band images of the galaxy 
sample. For nine out of the 20 galaxies, Sloan Digital Sky 
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Fig. 1. Example of the emission correction with GANDALF for the central spectrum of NGC 1358. Grey line shows 
the spectrum after being cleaned from emission lines, while the black spectra show the spectrum before. NGC 1358 is 
a galaxy from the sample with particularly strong emission lines, the fitting process applied to correct for the emission 
does a very good job even in this extreme case (Sect l4.11 for details). 



Survey! (SDSS) data was available. For NGC 1530 we use 
the data published in Zurita & Perez (2008), for NGC 2273 
the R — band image from Hameed & Devereux (1999) and 
for NGC 1169 the R - band from Knapen et al. (2004). 
Ellipse fitting to the light distribution is performed us- 
ing the ELLIPSE task within IRAF. As all the galaxies 
showed point-like nuclear regions, the center is fixed using 
the coordinates obtained by fitting a Gaussian to the nu- 
cleus. The position angle (PA) and the ellipticity (e) are 
left as free parameters in the fitting. The minimum ellip- 
ticity in the bar region has been adopted as the definition 
for the size of the bar. It has been shown (Michel-Dansac 
& Wozniak 2006) that this measurement is correlated with 
the corotation radius and, therefore, gives a bar structural 
size, it is also the method to calculate the bar size that is 
the least sensitive to dust absorption. Furthermore, it is a 
very simple criteria to reproduce and can be compared to 
other results in the literature. For a discussion about the 
differences in the bar size o b tained using different meth- 
ods see for example lErwin (2005). The bulge radius is 
considered as to the radius where the PA and ellipticity 
e starts showing the typical bar profile, i.e. increasing e 
and constant PA. Figure [2] and Fig. [3] show the results 
of the ellipse fitting for all the analyzed galaxies together 
with a R — band image, and an image of the residuals 
obtained subtracting the galaxy model to the R — band 



http:/ /www. sdss.org 



galaxy image. The radius of the bar, bulge and the ex- 
tent of our spectrocopic data are overplotted on the two 
figures. For the galaxies for which no SDSS data was avail- 
able, Fig. 2] shows 2MASS data together with the spectro- 
scopic data extent. Table [3] shows the derived bar size for 
the galaxies, with good imaging data from the SDSS and 
the literature, in comparison with the bar sizes obtained 
by other authors. There is good agreement between the 
values. The beginning of the bar coincides, in many cases, 
with a change on the trends of the line strength indices 
with radius (see Sec. 14. 4p . giving confidence in our size es- 
timations. The general morphology of the galaxies in the 
sample shows a large range of properties. Table Q] shows 
the detailed morphological types of all the galaxies. Some 
galaxies, independently of their type, have outer rings, 
some host strong bars and some are of intermediate AB 
type. We, therefore, cover a large range in general galaxy 
morphology. This allows us to test systematic behaviors 
associated to different morphologies. 

4.4. Derivation of stellar population parameters 

As we mentioned before (Section 14. 2p . to derive stellar 
population parameters we compare our line-strength in- 
dices with two different sets of stellar population models: 
TMB03 and V08. TMB03 use the fuel consumption theo- 
rem (Renzini & Buzzoni 1986) to evaluate the energetics 
of the post-main sequence phases. They use the stellar 
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Table 2. Log of observations. 



Table 3. Bar sizes. The bar size derivation is described 
in Section [4731 



Object Telescope Date 



PAjdeg) Exp. TimeG alaxy name 



NGC 1169 
NGC 1358 
NGC 1433 
NGC 1530 
NGC 1832 
NGC 2217 
NGC 2273 
NGC 2523 
NGC 2665 
NGC 2681 
NGC 2859 

NGC 2935 
NGC 2950 
NGC 2962 
NGC 3081 
NGC 4245 
NGC 4314 

NGC 4394 
NGC 4643 

NGC 5101 



INT 
INT 

2.3m,SSO 

INT 

INT 

2.3m, SSO 

INT 

INT 

2.3m,SSO 

INT 

INT 

2.3m,SSO 

INT 

INT 

2.3m,SSO 

INT 

INT 

INT 
INT 

2.3m,SSO 
2.3m, SSO 



09/02/2007 

10,11/02/2007 

03,06/02/2006 

08/02/2007 

12/02/2007 

6/02/2006 

07/02/2007 

09/02/2007 

01,02/02/2006 

02/02/2007 

08/02/2007 

12/02/2007 

5,6/02/2006 

11/02/2007 

10/02/2007 

01/02/2006 

07,08/02/2007 

09/02/2007 

12/02/2007 

11/02/2007 

10/02/2007 

1,5,6/02/2006 

3,5,6/02/2006 



95 

127 

96 

121 

171 

290 

116 

118 

54 

28 

162 

72 

135 

162 

175 

64 

135 

147 

237 

143 

134 

134 

122 



3h 

3h 

2.5h 

3h 

3h 

2h 

3h 

3h 

3.3h 

3h 

3h 

3h 

3h 

3h 

3h 

2.5h 

3h 

2.5h 

3h 

3h 

3h 

3.3h 

3h 



Bar size (arcsec) 
(other studies) 



Bar size (arcsec) 
(our study) 



tracks of Cassisi, Castellani & Castellani (1997), Bono et 
al. (1997) and Cassisi et al. (1999), and for the more metal 
rich models the ones by Salasnich et al. (2000). V08 uses 
the isochrones of Girardi et al. (2000) for all metallicities, 
and an empirical calibration to transform the models from 
the theoretical to the observational plane. By comparing 
with two different sets of models we take into account 
some of the uncertainties coming from the variety of in- 
gredients chosen by different authors. Below, we give a 
brief explanation of such possible uncertainties. 

V08 models provide synthetic spectra of single stellar 
populations with a given age, metallicity and solar chem- 
ical abundances ratios (the ones of the empirical library) . 
We know that, in some systems, the ratio between differ- 
ent elements do not match the solar values. For example, 
in giant elliptical galaxies, [Mg/Fe]>0.0 (e.g. O'Connell 
1976; Faber et al. 1992). It is important to take this effect 
into account because the presence of differences in the 
chemical abundance partitions affects the age measure- 
ments (e.g. TMB03). Several groups have characterized 
the variation of line-strength indices for single stellar pop- 
ulations to variations of individual chemical abundances 
(e.g. Tantalo et al. 1998; Trager et al. 2000; Proctor & 
Sansom 2002; TMB03). In order to do this, we follow the 
prescriptions given in Trager et al. (2000) with the cor- 
rections for negative indices given in Thomas et al. (2003) 
using the response functions by Korn et al. (2004). We 
make the assumption that all the a-elements are equally 
enhanced with the exception of Ca, which is assumed to 
follow Fe-peak element abundances. C is also assumed to 



NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



1169 

1358 d 

1433^ 

1530 e 

1832 / 

2217 i,;/;i 

2273 a 

2523 

2665 

2681 a ' 1;c;9 

2859 i,;c 

2935 /;i 

2950 a ' 3 

2962 a 

3081 b;9;c 

4245" 

4314 a ' 3 

4394^ 

4643 a 

5101 / 



25 

90, 66 

50 

17 

50, 37, 43 
14 



50, 25, 29 

34, 46 
25, 39 
24, 38 
29 

35, 41, 38 
37 

67, 75 
52, 35 
50 



2!) 



09 



21 



23 

48 

44 
45 
35 
59 
92 
56 
67 



50 



(1) N GC 2681 is a triple barred system. Erwin (2005) gives the 



value of the prim ary bar , whil e we give here the value of the 
secondary bar (a) lErwinl (|2005T l bar size given as the radius of 
ma ximum ellipti city in the bar region 

(b) lErwinl l|2004) bar size given at the radius of maximum el- 
lipt icity in the bar region 

(c) IWozniak et all (ll995T l. they give the projected bar size at 
the radius of minimum ell ipticity 

(d) iMulchaev et "ah! (|l997l ). they give the bar size as the radius 
at the minimum of elli pticity 

(el iRegan et all (Il996l) 

(f ) I Jungwiert et al. I" ifl997l ). bar size given as the radius of max- 
im um ellipticity in the bar region 

fgl lFriedli et~al] <|l996h . they give the projected bar size as the 
radi us at the ellipticity minim um 

(h) lGadotti fc de Souzal (|2006n . bar size given as the radius of 
m aximum elliptic i ty in the bar region 

(i) iMartin fc Rovl (l995l ) , bar size given by visual inspection of 
photographic plates 



be enhanced. Therefore, we define two groups: The en- 
hanced group (C, N, O, Ne, Na, Mg, Al, Si, S and Ar) and 
the depressed group (Fe, Ca, Cr and Ni). This matches the 
model 4 in Trager et al. (2000). We assume a linear rela- 
tionship between the [Z/H], [Fe/H] and [E/Fe] (where E 
represents the mass of all elements with enhanced abun- 
dance ratios) of the form [Z/H] = [Fe/H] + 0.929 [E/Fe] (see 
Trager et al. 2000) and all the variations of [E/Fe] are com- 
pensated by variations of the Fe-peak group in such a way 
that [Z/H] stay constant. 

The changes of the indices to different chemical com- 
position is calculated as (see Trager et al. 2000): 



s i = U( 1 + R ^( x ^ /H]/0 - 3 -^ 

-in 



(1) 
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where i?o.3(A"i) is the metallicity dependent response func- 
tion by Korn et al. (2004) to a variation of [X,/H] of the 
element i. 

The models of TMB03 already incorporate a similar 
treatment to take into account the differences in the chem- 
ical abundances ratios. We use the original version of the 
models with [C/Fe] and [N/Fe]=0 and [Ca/a]=0, but this 
choice does not affect our results. 

To obtain stellar population parameters we use 4 in- 
dices, namely; US a, H"fA, Fe4383 and Mgb and follow a 
multi-index approach as described in Proctor & Samson 
(2002). We avoid using H/3 because this is the Balmer in- 
dex most affected by the presence of emission. The errors 
are estimated by performing 100 Monte Carlo simulations 
where each index was perturbed by a gaussian distribution 
with a equal to the index error. 

The average deviations of the measured indices from 
the indices predicted by the model for the derived val- 
ues of age, metallicity and [E/Fe], are shown in Fig. IC.ll 
Despite the fact that we have only used 4 indices for the 
fit, it can be seen that most of the indices are repro- 
duced within the errors by the stellar population param- 
eters derived with this method. Some of the indices that 
are worse reproduced for some galaxies are the broad in- 
dices Mgi, Mg 2 , CNi and CN 2 . This is understandable 
as those indices are very sensitive to the flux calibration. 
As we have not corrected from internal extinction in the 
galaxy, the shape of the continuum in our spectra is af- 
fected by reddening. For the rest of the indices, the effect 
of dust is almost negligible (see, e.g. MacArthur 2005). 
In order to analyze the possible correlations between the 
results from the stellar population analysis and other pa- 
rameters, we fit a slope to the derived age, metallicity 
and [E/Fe] distributions with a simple linear regression 
weighted with the errors in the y-direction. Table [5] shows 
the fitted slope together with the errors. For four galaxies 
out of the 20, the obtained gradients are not reliable due 
to large fitting errors and these galaxies are discarded in 
the plots and the analysis regarding the gradient values 
(NGC 1530,NGC 3081,NGC 4314,NGC 2935), we include 
them in Tab. [5] 

To check the dependence of our results to the choice of 
stellar population model, we repeat the calculations using 
the new models by V08. Fig. [3] shows the comparison of 
the age, metallicity and [E/Fe] gradients (the slope of the 
fit) using different models. As can be seen, the comparison 
is very good in general. We have checked that none of our 
results is affected by our choice of model. 

An important concern about the obtained results is 
that a non-perfect emission correction could still be af- 
fecting our derived indices. The effect of a poor emission 
correction is to make the age older and, therefore, because 
the errors in age and metallicity are correlated, the metal- 
licity is lowered. It is difficult to estimate an error due to 
emission correction that takes into account the systematic 
effects as, for example, a template mismatch. Therefore, 
to test if our results are biased due to a poor emission sub- 
traction, we re-calculate the ages and metallicities using 8 



indices non-affected by emission: G4300, Fe4383, Ca4455, 
Fe4531, C4668, Mgl, Fe5270, Fe5335. Basically we use all 
the indices except for the ones that may be affected by 
emission. It has been shown that even without Balmer 
lines, it is possible to obtain a reliable age (although with 
larger errors) (e.g. Chillingarian et al. 2008). 

The differences between the mean gradients using 
these set of indices and using the previous set of 4 indices 
are shown in Fig.[6j As can be seen from this figure, there is 
very good agreement between the results obtained using a 
different set of index combination giving confidence in the 
derived results. However, for the [E/Fe] gradients there 
are three galaxies for which the values seem to slightly 
differ from the two methods. 

Before continuing with the presentation of the results, 
we want to remind the reader that the stellar popula- 
tion parameters derived here are single stellar population 
(SSP)-equivalent parameters and can be interpreted as 
mean values weighted with the light of individual stars. 
In the case of the age, the mean values are going to be 
strongly biased towards the age of the youngest compo- 
nents. 

5. Results 

5.1. Kinematics 

Figure [8] shows the line-of-sight position- velocity diagrams 
and velocity dispersion along the bar major-axis for the 
galaxies in our sample. Close examination of these dia- 
grams can help to reveal the nature of the structures in 
the bar region. All the galaxies present disk like structures 
from their kinematics. These structures are characterized 
by a central velocity dispersion minima (er drops) , central 
a plateaus and/or a clear rotating inner disks or ring as 
derived from the line-of-sight position-velocity diagrams. 
We define the central a as the maximum velocity disper- 
sion in the bulge region. We have chosen this dispersion 
value to avoid being dominated in some galaxies by the 
low central velocity dispersion that could be indicative of 
nuclear star formation. Therefore, the maximum velocity 
dispersion should be closer to the velocity dispersion of 
the bulge. 

Most of our galaxies show " double- hump" rotation 
curve (Bureau & Athanassoula 2005; Chung & Bureau 
2004), that is, the rotation curve first rises rapidly and 
reaches a local maximum, then drops slightly and then 
raises again slowly. These profiles were detected in N- 
body simulations of edge-on barred disks by Bureau & 
Athanassoula (2005). They appear in strong bars and are 
detected in small viewing angles. The fact that these fea- 
tures are detected in pure N-body simulations indicate 
that they could arise naturally from the orbital structure 
of a barred disk and therefore gas dissipation is not in- 
dispensable. Central velocity dispersion minima occur in 
70% of our sample. Galaxies not showing a cr-drop are : 
NGC 2859, NGC 2935, NGC 2950, NGC 2681, NGC 2217, 
NGC 1169 and NGC 1358. Although, the NGC 2950, 
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NGC 2217 and NGC 2859 shallow and wide central dip 
could be considered as well as a cr-drop. Chung & Bureau 
(2004) found c-drops in 40% of their barred sample, while 
Peletier et al. (2007) found it in 50% of their sample of 24 
barred and unbarred galaxies. Marquez et al. (2003) re- 
ported that 12 out of 14 of the analyzed galaxies with 
CT-drops showed an inner disk from their HST images. 
The standard explanation for these drops is the presence 
of cold central stellar disks originating from gas inflow 
(Emsellem et al. 2001; Wozniak et al. 2003). Wozniak & 
Champavert (2006) predicted from dynamical simulations 
of drops in barred galaxies that these features can last for 
more than 1 Gyr if the central region is continuously fed 
by fresh gas leading to a continuos star formation activity. 
In their simulations, continuous star formation is needed 
to replace the heated particles by new low-cr ones. They 
predict that as the star formation disappears, the drop 
will disappear as well. However, this disappearance is not 
instantaneus, they find that a cr-drop dissapears at a lin- 
ear rate of 10 kms -1 in 500 Myr, when star formation is 
switched off. On the contrary, Athanassoula & Misiriotis 
(2002) found these features in N-body dissipationless sim- 
ulations, and Bureau & Athanassoula (2005) showed that 
they can arise from the orbital structure of strongly barred 
galaxies. These authors noted that, in their simulations of 
strong barred galaxies, the width of the central flat section 
of the velocity dispersion profiles was roughly the same as 
that of the rapidly rising part of the rotation curve, and 
suggest a common origin (xl orbits). However, these N- 
body simulations have no stellar populations predictions 
and, therefore, it is difficult to compare to our observa- 
tions. These two possible cr-drop origins (cold dissipative 
disk vs. bar orbital stellar structure) can be separated by 
studying the age in the central parts of the galaxies with 
er-drops. If the cr-drop has to be maintained by continuous 
gas replenishment and star formation, we would expect a 
very young population dominating the central regions. On 
the other hand, if the drops are due to an old stellar struc- 
ture, we would expect the central parts to be dominated 
by an old stellar component. We will present the stellar 
population of the central parts of these galaxies in the 
second series of the paper, but we can already advance 
that, while most of these galaxies have, indeed, young 
components and star formation in their centers, eight 
galaxies (NGC 1358, NGC 2523, NGC 2859, NGC 2962, 
NGC 4245, NGC 4314, NGC 4643 and NGC 5101) have 
old central luminosity weighted ages (> 3 Gyr) despite 
having central dips or plateaus. However, this fact does 
not, necessarily, imply that the ages of the cold compo- 
nent are old, since we are measuring luminosity weighted 
ages and a bulge component could be contributing to the 
derived old age. It would be necessary to disentangle the 
populations coming from both components. This will be 
done in detail in a forthcoming paper addressing the cen- 
tral regions of barred galaxies. 

Eight out of the 20 galaxies are classified as double 
barred (Erwin 2004) NGC 1433, NGC 2217, NGC 2681, 
NGC 2859, NGC 2950, NGC 2962, NGC 3081, NGC 4314. 



They do not show any difference in their kinematic pro- 
files, although all of them show a plateau in their central 
stellar velocity dispersion. 

5.2. Mean stellar population parameters in the bar 
region 

In order to relate the bar properties with their stellar pop- 
ulation parameters, we first calculate the mean stellar age, 
[E /Fe] , and metallicity of the bar by averaging the derived 
values along the bar and weighting with the errors in these 
parameters. At this point, we would like to clarify the dif- 
ference between what we mean by the age of the bar and 
the age of the bar stellar population, the later is the mean 
age of the stars as derived from the stellar populations, 
and the age of the bar is, for us, the age constrains that 
we obtain from the stellar population analysis, mean val- 
ues and gradients of the ages and metallicities, to the time 
of bar formation. The mean values are reported in TableSJ 
The first thing that can be seen is that our sample of bars 
exhibits a large dispersion in the mean ages of their stel- 
lar populations. The differences in the mean ages in the 
bars between the different galaxies can be as large as 10 
Gyr, with some galaxies showing a mean bar age larger 
than 10 Gyr and others younger than <1 Gyr. Figured] 
shows the relationship between the mean stellar age, the 
mean Z and the average [E/Fe] in the region of the bar 
against the maximum central velocity dispersion. As can 
be seen, there is a trend in which galaxies with smaller cen- 
tral velocity dispersion are also the ones showing younger 
mean- ages. In the same way, galaxies with lower maxi- 
mum central velocity dispersion tend to have lower values 
of [E/Fe] and metallicity. In this Figure, we have separated 
the galaxies hosting a nuclear AGN from those without, 
there is no significant difference in the distribution of the 
two samples. 

This is the first time that stellar population parameters 
are derived for the bar region in galaxies. It is interesting 
to notice that the trends shown in Fig. [9] are similar to 
those found for galaxy bulges (Proctor & Sansom 2002, 
Morthy & Holtzman 2006) with values of the metallicity 
and [E/Fe] slightly lower than those found in the bulges. 
This result points out to an intimate link between the 
evolution of the bar and the bulge. 

Gadotti & de Souza (2006) also found, using color as 
a proxy for ages, a large difference in the ages of the stars 
in a sample of 18 barred galaxies. These authors found 
a relation between the color of the bar and the morpho- 
logical type of spiral, being more likely that a later type 
spiral host a younger bar. We did not find any correlation 
between the mean age of the stars in the bars and the 
morphological type of the galaxies, although we do not 
have late spirals in our sample. 

These authors also found a weak correlation between 
the mean age of the bar and the length of the bar, with 
older bars being also longer. We did not find any correla- 
tion between these two parameters. As a cautionary note, 
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Table 4. Mean stellar population parameters. Mean ages 
(2), mean metallicities (3) and [E/Fe] (4) (see Section l4~4l 
for an explanation on the parameter derivation). 



Table 5. Metallicity and age gradients 



Galaxy name 


< lo 


I age > 


<[Z/H]> 


<[E/Fe]> 


NGC 


1169 


1 


172 


± 





038 


0.069 ± 0.032 


0.112 ±0.029 


NGC 


1358 


1 


043 


± 





048 


0.000 ± 0.034 


0.225 ± 0.029 


NGC 


1433 





266 


± 





031 


0.005 ± 0.041 


0.138 ±0.038 


NGC 


1530 





078 


± 





027 


-0.023 ± 0.050 


0.267 ±0.073 


NGC 


1832 





153 


± 





015 


-0.184 ±0.028 


0.002 ± 0.024 


NGC 


2217 





344 


± 





045 


0.644 ± 0.046 


0.336 ± 0.030 


NGC 


2273 





241 


± 





020 


0.064 ± 0.030 


0.111 ±0.025 


NGC 


2523 





530 


± 





062 


0.313 ±0.054 


0.136 ±0.031 


NGC 


2665 





334 


± 





076 


-0.606 ± 0.092 


-0.010 ±0.076 


NGC 


2681 





439 


± 





029 


-0.602 ±0.039 


0.015 ±0.036 


NGC 


2859 


1 


022 


± 





029 


-0.064 ±0.021 


0.238 ±0.021 


NGC 


2935 





374 


± 





056 


-0.192 ±0.059 


0.114 ±0.052 


NGC 


2950 


1 


053 


± 





025 


-0.311 ±0.019 


0.114 ±0.052 


NGC 


2962 





920 


± 





035 


0.029 ± 0.026 


0.135 ±0.022 


NGC 


3081 





950 


± 





157 


-0.225 ±0.124 


-0.060 ±0.096 


NGC 


4245 





774 


± 





045 


-0.095 ±0.031 


0.085 ± 0.028 


NGC 


4314 





297 


± 





031 


0.066 ± 0.038 


0.105 ±0.032 


NGC 


4394 





956 


± 





033 


-0.202 ± 0.024 


0.131 ±0.025 


NGC 


4643 





741 


± 





013 


0.321 ± 0.009 


0.126 ±0.006 


NGC 


5101 





477 


± 





021 


0.182 ±0.213 


0.113 ±0.015 



Galaxy name 


dfZ/Hl /dr 




NGC 


1169 


n noo _i_ n nno 

U.UZrO ZTZ U.UUO 


n noo4 4- o ooors 

U.UUU1 ZTZ U.UUUU 


NGC 


1358 


U.UitJ ZTZ U.U-L-L 


— 009Q 4- 09^^ 

U.UU^jC? ZTZ U.U^UO 


NGC 


1433 


n nnzi -I- n nn^i 

U.UU4: ZTZ U.UUl) 


u.u±4:( zrz u.uuyo 


NGC 


1530 


a nco ±nmo 

U.UtJ^r ZTZ U.U-LZr 


— 01 4 1 } 4- 0070 
u.ui^o zrz u.uuiu 


NGC 


1832 


_n nnc _i_ n nno 
u.uuu zrz u.uuo 


— 0989 4- 01 1 

U.U^O^ ZTZ U.U-L-LU 


NGC 


2217 


n noo _i_ n m o 
L/.vjZrO zrz 


— 0^94 4- 01 40 
\j .\jo zrz u.uj-tcU 


NGC 


2273 


o ni q -\- o ni i 

U.yJ-LcJ ZTZ U.UXJ. 


— 01 99 4- OOfi^ 
zrz u.uuuo 


NGC 


2523 


_n noi _i_ n p.97 


0000 4- 09QQ 


NGC 


2665 


_n 087 + 0014 

U.UOI ZTZ U.U-L4: 


040^ 4- 01 ^0 
u.uiuu zrz u.uiiju 


NGC 


2681 


— 0Q7 4- o^ 

\J.VJij i ZTZ U.UOO 


C\C\RR + 01 so 
u.uuoo zrz u.uiou 


NGC 


2859 


— n nns 4- n nnfi 

U.UUO ZTZ U.UUU 


007Q 4- 0089 

U.UU ( a ZTZ U.yJyJO^ 


NGC 


2935 


n of. a + n run 

U.UU4: ZTZ U.U4:U 


_n i 071 4. n n^OO 
u.±o( 1 zrz u.uouu 


NGC 


2950 


_n ooo _i_ n 00A 

U.UUU ZTZ U.UU4: 


0047 4- OO^Q 
u.uui/ zrz u.uuoy 






U.U^-L ZIZ U.U-LU 


— 09^7 4- 01 ^8 
u.u^u ( zrz u.uioo 


NGC 


3081 


-0.587 ±0.556 


-0.0661 ± 0.6848 


NGC 


4245 


-0.010 ±0.005 


0.0337 ± 0.0098 


NGC 


4314 


-0.014 ±0.005 


0.0184 ± 0.0096 


NGC 


4394 


0.023 ± 0.005 


-0.0246 ± 0.0066 


NGC 


4643 


-0.004 ±0.002 


0.0063 ± 0.0044 


NGC 


5101 


0.008 ±0.002 


-0.0268 ±0.0012 



it is not a good assumption to consider broad-band colors 
as a probe for age, assuming a constant metallicity for all 
the galaxies, as we have shown here that bars in different 
galaxies show a wide metallicity distribution. s 

\ 

M 

5.3. Stellar populations along the bar 



As shown in Sec. 14.41 we fit a slope to the derived age, 
metallicity and [E/Fe] distributions with a simple lin- 
ear regression weighted with the errors in the y-direction 
(Table[5]). We find that most galaxies show a gradient dif- 
ferent from zero in both age and metallicity. Fig. [7] shows 
the indices at the beginning and end of the bar fSec l4.3|l 
together with the predictions of the models by V08. The 
galaxies with the lowest mean metallicity values have, sys- 
tematically, younger and more metal poor populations at 
the end of the bar (these are the data points on the left- 
hand side of the grid in Fig. [7]). The more metal rich bars 
(data points on the right-hand side of the grid in Fig. [7|) 
show more complex behavior with most of the galaxies 
showing younger and more metal rich populations and the 
end of the bar. 

Figure [TO] shows the age and metallicity distribution 
profiles for the galaxies in our sample. We have not in- 
cluded in this analysis NGC 1530, NGC 3081, NGC 4314 
and NGC 2935 due to the large fitting errors. Seven galax- 
ies have metallicity gradients significantly (more than 3-cr 
significance) different from zero: NGC 1169, NGC 2217, 
NGC 4394, and NGC 5101 (positive) NGC 2665, 
NGC 2681, NGC 4245 (negative). 

We therefore, confirm the results from Perez et al. 
(2007), where we found, in a much more reduced sample, 



" ,-Vt 



200 
cj(km/s) 



Fig. 11. Metallicity difference between the inner and the 
outer bar region as a function of maximum central velocity 
dispersion. Notice that galaxies with a < 170 km s -1 tend 
to have negative gradient values. 



that some of the bars had positive metallicity gradients 
and negative age gradients. These negative age gradients 
are expected since, as it has been shown in numerical sim- 
ulations, there is an accumulation of young stellar popu- 
lation trapped on elliptical-like orbits along the bar, near 
the ultra- harmonic resonance (Wozniak 2007). Although 
expected, this negative age gradient should be further in- 
vestigated by the models to check whether the gradient 
can stay for longer than 3 Gyr, time at which the simula- 
tions where stopped and shorter than the age of the last 
star formation bursts found for our galaxies. 
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What we found more striking is the presence of posi- 
tive metallicity gradient in four galaxies, i.e. the younger 
stars at the end of the bar are also more metal rich. These 
gradients are difficult to understand in the context of our 
current understanding of bar formation (see Section 6) , if 
they are a consequence of the star formation ocurred dur- 
ing its formation, this would imply that there has been 
chemical enrichment and, therefore, that the star forma- 
tion lasted for a long period of time or that it has been 
very efficient. 

The galaxies showing a negative metallicity gradient 
(i.e. less metal rich in the outer parts of the bar) might 
be reflecting the original disk gradient. Only for a few 
galaxies the disk radial stellar gradients have been de- 
rived. Molla, Hardy & Beauchamp (1999) derived shal- 
low negative gradients for the Mg2 and Fe5270 indices 
in the discs of NGC 4303 (SABbc), NGC 4321 (SABbc) 
and NGC 4535 (SABc). Ryder, Fenner & Gibson (2005) 
obtained the radial disk profiles of the Mg2 and Fe5270 
indices for a sample of 8 late-type galaxies with similar 
results. In a recent work, Yoachim & Dalcanton (2008) 
have derived luminosity weighted ages an metallicities in 
the disk region of nine edge-on galaxies, finding radial neg- 
ative age gradients and almost no metallicity gradient. In 
order to follow up this issue, we have collected deep spec- 
troscopic information in the disk of our sample (Sanchez- 
Blazquez et al. 2009, in preparation). 

Figure QT] shows the relation between the metallic- 
ity gradient along the bar and the velocity dispersion. It 
can be seen that the galaxies with velocity dispersions 
higher than 170 kms" 1 tend to have positive metallic- 
ity gradients, while the lower velocity dispersion galax- 
ies show a larger dispersion in their metallicity gradient 
value. Similar behavior is found for the age gradient, where 
galaxies with large velocity dispersion tend to have a shal- 
low positive age gradient while for the lower velocity dis- 
persion galaxies the absolute value of the age gradient is 
larger, both negative and positive gradients. This seems 
to imply that the stellar population gradient is, somehow, 
related to the mass of the galaxy's bulge. To test whether 
there is a link between the mass of the bulge and the total 
galaxy mass, we have investigated the correlation between 
the bulge dispersion velocity and the total ii-band lumi- 
nosity (as a tracer of the stellar galaxy mass). There is, in- 
deed, a correlation between the maximum central velocity 
dispersion and the total iJ-band luminosity, which might 
link the observed metallicity gradients to the stellar mass 
of the galaxy. Further investigation needs to be done to 
address this issue. The galaxies with negative metallicity 
gradients have systematically, positive age gradients and 
have average bar-age values of the order of 2 Gyr (except 
for NGC 4245 with an average of 5 Gyr). 

5.4. Correlation with bar/galaxy parameters 

We want to check possible connexions between observa- 
tional structural bar parameters, that are thought to be 



important in the bar formation and evolution, with the 
stellar population parameters derived in this work. The 
bar-strength is determinant in the orbital structure of 
bars and it is thought to evolve as the bar ages (e.g. 
Athanassoula 2005). For this purpose, we have used the 
available bar-strengths obtained as the ratio of the tan- 
gential force to the mean radial force derived from the 
NIR images (Combes & Sanders 1981, Block et al. 2001, 

2004) . For 13 out of the 20 galaxies there are values for the 
bar-strength from the literature (Table QJ. No clear rela- 
tion is found between the bar-strength and the metallicity 
gradients. 

It has been suggested that older bars should get longer 
and thinner in time as their pattern speed decreases, in 
order to check if there is any systematic behavior be- 
tween the bar sizes and the metallicity gradients, we have 
checked the absolute bar sizes and the relative sizes of the 
bar with respect to their disk for the galaxies in our sam- 
ple. We find, in agreement with previous studies (Erwin, 

2005) , that the bar sizes are amongst the typical sizes of 
SO-Sab (between 1 and 8 kpc in our sample) and that 
the sizes correlate well with the disk size (Perez et al. 
2005, Erwin 2005). There is no significant trend between 
the size of the bar and the metallicity and age gradient. 
However, there might be a hint of correlation pointing to 
larger bars having positive age gradients, negative metal- 
licity gradients and an average young population. But, due 
to low number statistics and the large spread in the values, 
nothing conclusive can be said to this respect. However, 
it would be interesting to follow-up this point calculat- 
ing the bar-size for all the galaxies with new broad-band 
data to test if there is a significant trend of the bar stellar 
population parameters with bar-size. Rings are located 
near resonances induced by the bar component, from nu- 
merical simulations outer rings form in time-scales longer 
than 1 Gyr and therefore, they have been associated to 
and relatively aged system. However, in some simulations, 
rings appear very quickly, about 0.2 Gyr after the bar 
formation (Rautiainen & Salo, 2000). We do not see any 
correlation between the presence of an outer ring and the 
bar age. Some of the bars with young populations and 
negative metallicity gradients present outer rings, there- 
fore, the presence or absence of an outer ring is not a good 
discriminator of the bar age. 

There seems to be no clear correlation with morpho- 
logical type, although we do not cover a large range in 
morphological types (SO-Sbc), the galaxies with positive 
metallicity gradients in the bar region are not particularly 
different in morphology to the rest of the sample. 

5.4.1. AGN connection 

Since the bar is an efficient mechanism to drive material to 
the central parts of the galaxy it has been very much dis- 
cussed its role in feeding the central massive black hole and 
the bar-AGN connection. Recent statistical studies (e.g. 
Knapen, 2004) have concluded that there is not a statisti- 
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cally significant relation between the presence of an AGN 
in barred galaxies. Part of the problem in the lack of cor- 
relation between the AGN and the bar presence could be 
the different time-scales of the different processes. There 
have been recent claims, using broad-band photometry 
(Gadotti & Souza, 2004) that active galactic nuclei tend 
to have young bars, while galaxies with old bars do not 
generally host AGNs. These studies are based on the as- 
sumption that the change in broad-band galaxy color is 
due to age and not to metallicity as well. Because spec- 
tra provide a better way to discern between both age and 
metallicity, and since the indices are not affected by dust, 
we have carried out this project also to see whether there 
is any correlation between the bar ages and metallicities 
with the presence of an AGN. In Fig.[9]we have separated, 
with different symbols, galaxies hosting an AGN. We do 
not see any significant correlation with mean bar age and 
metallicity nor with the derived gradients. It might only be 
reflecting the small number statistics and a larger sample 
would be needed to properly cover the parameter space. 
Therefore, we do not reproduce the results of Gadotti & 
de Souza (2006). However, the correlation that these au- 
thors found is weak and they were using colors as a proxy 
for age. 

6. Discussion and Summary 

We have derived the kinematics and the distribution of 
ages and metallicities along the bar in a sample of 20 
barred galaxies to study the question whether bars are 
long-lived or weather they are destroyed and reformed in 
a few Gyrs. We also want to explore the relation between 
the age, and stellar populations of the bar and other global 
physical properties of the host galaxies. 

From the study of the kinematics along the bar, we 
have found that all the galaxies in the sample show a 
disk-like component in their centers, showing as dips or 
plateaus in their central stellar velocity dispersion and as 
stellar rotating disks. Therefore, this implies that bars, as 
expected, have a strong influence in the building up and 
later evolution of the central component. This aspect will 
be studied in more detail in Sanchez-Blazquez, Perez and 
Zurita (2009, in prep) where we compare the bulges ages 
and metallicities of barred and unbarred bulges of other- 
wise similar galaxies. 

We have derived ages and metallicities from Lick/IDS 
line-strength indices of single stellar population (SSP). 
It is clear that the true star formation history in these 
objects will be more complicated than a simple SSP. 
However, a lot of information can be derived from the 
SSP analysis. The SSP-equivalent age of a composite stel- 
lar population (CSP) is dominated primarily by the age 
of the youngest component and the mass fraction of the 
different populations. The SSP-equivalent metallicity re- 
flects a V — band luminosity weighted chemical composi- 
tion (Serra and Trager 2007). 

As bars form from stars already in the disc, we have to 
distinguish between stellar population gradients that may 



have already been present when the bar formed and those 
that formed during the bar formation or afterwards. There 
are not many studies analysing the stellar population of 
discs outside the MW. In particular, for early-type spi- 
rals, only broadband colours have been traditionally used 
to estimate the ages and metallicities of the stellar discs. 
However, stellar parameters are very difficult to obtain 
with broad band colours. Recently, Yoachim & Dalcanton 
(2008) have measured Lick indices in the thin and thick 
disc of edge-on galaxies. They found a very strong age 
gradient and a null or negative metallicity gradient de- 
pendending on the index they were using. In the MW a 
negative age and metallicity gradient are measured. We 
will assume, for our discussion, that this is the general 
trend on the discs before the formation of the bar. 

We have found three different types of bars according 
to their metallicity and age distribution along the bar: 

1) Bar with negative metallicity gradient. These bars 
show a mean young/intermediate population (< 2Gyr) 
and have amongst the lowest central stellar velocity dis- 
persion of the sample. They also tend to have positive age 
gradients (i.e. younger populations at the beginning of the 
bar). In fact, all the galaxies showing these trends posses 
circumnuclear regions with young population ( 1.5 Gyr). 
The negative gradient can be understood as the original 
gradient of the disc. This is because a small fraction of 
young stars can bias the mean-age measured with SSP to- 
wards very small values. However, the mean metallicity 
is more sensitive to the more massive component. In this 
picture, the circumnuclear region would be strongly af- 
fecting the mean age at the beggining of the bar (lowering 
it) but its influence on the mean metallicity would not be 
so strong. Since one would expect the original metallicity 
gradients to flatten by the velocity dispersion of the bar 
(Friedli et al. 1994), the presence of metallicity gradients 
may be indicating that they formed recently. 

2) Bars with null metallicity gradients. The galaxies 
that do not show any gradient in their metallicity distri- 
bution along the bar have negative age gradients, that is, 
their ends of the bar are younger. These galaxies present 
gradients in both metallicity and age compatible with the 
models that studied mixing of metals and the distribution 
of stellar population ages in bars using numerical simula- 
tions (e.g Wozniak 2007; Friedli et al. 1994). They predict 
low age regions at the end of bars due to the accumula- 
tion of young stellar populations trapped in orbits near 
the ultra-harmonic resonance. During the first Gyr, in the 
simulations, most of the star formation happens. We ob- 
serve that, although, the population at the end of the bar 
is younger than the average bar value, it is still larger than 
3 Gyr. This fact might be indicating that the bar in these 
galaxies is at least 3 Gyrs old. 

3) Bars with positive metallicity gradients. These 
galaxies are the ones with older mean ages in their stel- 
lar populations. They are also the ones with more mas- 
sive bulges (larger a). In principle, this positive metallic- 
ity gradients are difficult to understand in light of our 
current knowledge of formation and evolution of bars. 
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Hydrodynamical simulations indicate that, even if metal 
rich population can be formed at the end of the bar, the 
diffusion timescale in this region is very short, of the or- 
der of 50 Myr (see figure 3 in Wozniak 2007) . The metal- 
licity gradient in disc galaxies is believed to be negative 
(see references above) and, therefore, it is not likely that 
this gradient is a relic of the original metallicity gradi- 
ent on the disc. We do not have a definitive answer to 
the origin of this metallicity gradients. One posibility is 
the presence of an inner ring at the beggining of the bar 
with enriched SF that could be populating the extremes 
of the bar in the same fashion as we explained for the 
presence of young stellar populations in the outer parts of 
the bar (case 2). However, we do not see this effect in all 
the galaxies hosting a nuclear young ring. Furthermore, 
the luminosity weighted mean age of the stars at the beg- 
gining of the bar is older than that in the external parts 
which is not expected if the bar ends are being populated 
with an enriched population coming from the inner ring. 
Projection effect of the old metal-rich population of the 
bulge on the bar region is not a likely scenario, since, at 
it can see from the kinematics, the outer part of the bars 
are dominated by a cooler stellar component. 

It has been claimed that the star formation in bars 
follows a sequential pattern, where younger bars (< lGyr) 
would show star formation along the entire bar and later 
the star formation would be concentrated on the bar-ends 
and the nuclear region (Phillips 1993; Martin & Friedli 
1997; Martin & Roy 1995; Friedli & Benz 1995; Verley 
et al. 2007). We have checked the galaxies for emission 
lines coming from HII regions to see if the ones showing 
negative metallicity gradients present systematically HII 
along their bars. 

NGC 2665 and NGC 1530 show nebular emission along 
their bars belonging to HII regions, which reflects very re- 
cent star formation events. Unfortunately, we do not have 
high S/N information on the stellar indices for the whole 
bar in NGC 1530 that would allow us to establish the 
stellar population age and metallicity. NGC 2665, on the 
other hand is one of the galaxies that show a negative 
metallicity and a positive age gradient along its bar. This 
is compatible with the idea of this sequential star forma- 
tion pattern if we assume that these bar have formed less 
than 3 Gy ago. However, other bars such as, NGC 4245 
and NGC 2681, although belonging to the same case as 
NGC 2665, do not present any HII regions along their 
bars. 

The lowest average bar value of the a/Fe, here ex- 
pressed as [E/Fe] (Trager et al. 2000), tends to have lower 
central velocity dispersion values (Fig. [9|), these galaxies 
have also in general a lower average age (Fig. [9|). The 
galaxies with higher central velocity dispersion stopped 
forming star earlier than galaxies with lower velocity dis- 
persion. However, this is not an indication that the galax- 
ies with lower central velocity dispersion host a younger 
(i.e. a more recently formed) bar. 



7. Conclusions 

We have measured, for the first time, line-strength indices 
along the bar of a sample of 20 early-type spiral galaxies 
and SO. We find that the mean values of age, metallicity 
and [E/Fe] correlate with central velocity dispersion in 
a similar way to that of bulges, pointing to a intimate 
evolution of both components. Galaxies with high stellar 
velocity dispersions (>170 kms -1 ) host bars with old stars 
while galaxies with lower central velocity dispersion show 
stars with a large dispersion in their ages. 

We find, for the first time, gradients in age and metal- 
licity along the bars of all the galaxies. We have found 
three types of bars according to their metallicity and age 
distribution along the radius. Most galaxies could be com- 
patible with having formed more than 3 Gyr ago, although 
a few of them show characteristics compatible with having 
been formed less than < 2 Gyr ago. These derived gradi- 
ents place strong constrains on models of bar evolution. 

We have explored the possible correlation of the bar 
age with other bar/galaxy parameters, such as bar- 
strength, bar size, morphological features, and AGN pres- 
ence without finding any significant trend. All the galaxies 
show disk-like central components, implying a strong role 
played by bars in the bulge secular evolution. Among these 
disk-like central components we find some galaxies hosting 
an old, luminosity weighted, stellar component (> 3 Gyr) 
and others with a recent (< 1 Gyr) burst. Therefore, it 
would interesting to separate the bulge an inner disk pop- 
ulations to analyse their stellar population parameters in 
order to investigate the origin of the u-drops. This will be 
presented in a forthcoming paper. 
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Appendix A: Line-strength distribution 

Fig. I A. II shows the line-strength gradients in the bar region 
for all the galaxies. 
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Fig. 2. Morphological analysis of the sample galaxies from SDSS images. For every galaxy it is shown the ellipticity 
distribution, position angle distribution the i?-band image, the red circle shows the extent of the data and the black 
circle the position of the bar end, and the residual image. 
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Fig. 3. Continuation of Fig [2] 
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Fig. 4. 2MASS images of the galaxies without SDSS imaging. The red circles show the extent of the spectroscopic 
data 
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grad (Age) TMB03 grad[Z/H] TMB03 grad[E/Fe] TMB03 

Fig. 5. Comparison between the mean gradients measured using TMB03 and V08 models. We have carried out a 
comparison of the metallicity and age gradients using two different models to asses the dependance of the results with 
the model choice. As it can be seen, the comparison is very good in general. 
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Fig. 6. Comparison between the age, metallicity gradients obtained using 4 indices or a combination of 8 indices not 
affected by emission (see text for details). There is very good agreement between the results obtained using different 
sets of index combination. 




1 1.5 2 2.5 3 3.5 4 2 4 

[MgbFe] Fe4383 



Fig. 7. [MgFe] and Fe4383 indices vs. US a measured in the inner (filled symbols) and outer (open symbols) parts of the 
bar. Overplotted are the models by Vazdekis. Models are shown with ages = 1,1.41,2.00,2.82,3.98,5.62,7.98,11.22,15.85 
(Gyr) from top to bottom, and metallicities, [Z/H] = -1.68,-1.28,-0.68,-0.38,0.0,0.2 from left to right. 
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Fig. 8. Position- velocity diagrams and velocity dispersion profiles along the bar major-axis for our sample of galaxies. 
Dashed lines indicate the beginning of the region dominated by the bar. 
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Fig. 8. Position- velocity diagrams and velocity dispersion profiles along the bar major-axis for our sample of galaxies. 
Dashed lines indicate the beginning of the region dominated by the bar 
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Fig. 9. Mean age, metallicity and [E/Fe] in the bar region against maximum central velocity dispersion. The diamonds 
correspond to the galaxies classified as hosting and AGN, the circles are the galaxies with non-active nuclei. Notice 
that galaxies with lower maximum central velocity dispersions are also the ones with younger ages and lower Z and 
[E/Fe] values, these trends are similar to those found for the bulge region 
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Fig. 10. SSP-equivalent age and metallicity along the radius. Dashed lines indicate the beginning and the end of the 
bar region. A linear fit to the bar region is also plotted. 
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Fig. 10. SSP-equivalent age and metallicity along the radius. Dashed lines indicate the beginning and the end of the 
bar region. A linear fit to the bar region is also plotted. 
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Fig. 10. SSP-equivalent age and metallicity along the radius. Dashed lines indicate the beginning and the end of the 
bar region. A linear fit to the bar region is also plotted. 
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Fig. 10. SSP-equivalent age and metallicity along the radius. Dashed lines indicate the beginning and the end of the 
bar region. A linear fit to the bar region is also plotted. 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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index 



Offset (MILES) Offset (Lick) 



this work-miles this work-Lick 



H5 A 

m F 

CNi 
CN 2 



-0.18 0.00 
0.000 0.15 



0.000 -0.011 



0.000 0.00 

0.0 0.00 

0.2 0.0 

0.0 0.00 

0.0 0.00 

0.21 0.00 

0.09 -0.25 

0.12 0.00 



Ca4227 



G4300 
F4383 



Ca4455 
Fe4531 



C4668 
H/3 



-0.75 -0.50 
0.00 0.163 
0.00 0.00 



Fe5015 



Mgl 
Mg2 
Mgb 



-0.017 -0.030 

0.0 -0.030 

0.0 0.00 

0.0 0.0 

0.0 0.0 

0.0 0.0 



Fe5270 
Fe5335 
Fe5406 



Table B.l. Mean offsets measured in the stars in common 
between our observing runs and MILES (first column) and 
Lick/IDS (second column). 



Appendix B: Calibration to the models 
spectrophotometric systems 

The original Lick/IDS spectra were not flux-calibrated by 
means of a flux-standard stars but normalised to a quartz- 
iodide tungsten lamp. The resulting continuum shape can- 
not be reproduced and causes small offsets for the in- 
dices with a broad wavelength range. Furthermore, the 
wavelength dependent resolution of the Lick spectrograph 
is also difficult to reproduce, which cuases differences in 
other indices with larger dependence to resolution. In or- 
der to compare our measured indices with models based 
on this library we have to calculate these offsets. To cal- 
culate them, we observed 11 and 35 stars in common with 
the Lick/IDS library. By comparing the indices measured 
in our stars with those in the Lick/IDS database for the 
same objects, we derived mean offsets for all the indices 
observed. Fig. IB. II shows this comparison. We did not find 
any systematic difference between the offset obtained with 
the stars of the first and second run and , therefore, we 
calculated just one offset for the two runs. The final offsets 
are listed in Table [B~T1 

Appendix C: Stellar Population fits 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. A.l. Line-strength distribution in the bar region for all the galaxies 
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Fig. B.l. Comparison of the indices measured in the stars in common between this work and the Lick/IDS library. 
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Fig. C.l. Deviations from best fit values for data compared to the models of TMB03. The indices that have been used 
for the fit are represented with filled symbols. 



